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SECTION SEVEN - NUMERICAL 
. M-ETHODSO'F SOLUTION 

USE OF NUMERICAL ANALYSIS IN INTERNAL BALLISTICS 
^Various variable quantities possessing definite physical 
significance usually take part in processes which occur in nature 
or are considered in technology. In this connection, numerical 
variations of one or more quantities are accompanied by or 
.associated with variations of other variables. Thus there always 
\exists a definite functional relation between the variable quantities 
vnder consideration. This functional relation may be expressed 
by means of three methods, such as tabulations, diagrams, and 
^formulas. In .the vast majority of procesces, especially those 
encountered in technology, this relation is expressed by the aid 
ofjtables or diagrams obtained directly from experiment or from 
observation of the process, whereas the formulas appear only after 
subsequent analysis of the results obtained, and then only in the 
case of the simplest processes. It is thus apparent that the most 
^natural means of expressing a functional relation between variable 
quantities representing in their totality the process under 
investigation is a tabulation; this is especially true when such 
r a process is being utilized directly for technological purposes, 

I_i« which case a formula or even a diagram will not serve the purpose, 
:and only numerical values of the variable quantity considered to 
be of. primary importance on the basis of practical considerations 
must fee had. 

-Numerical analysis must thus serve as a means for studying and 

making practical utilization of the functional relations between the 
F-7S-7327-RE 612 


CHARLES A. MEYER & CO. Inc. 













variable quantities involved. J 

Ordinary differential equations may bo approximately integrated j| 

by means of any of the known methods of which there are a great many. ,| 

s 

These methods include the following*. ".3— < 

1) Expansion in a Taylor's series in powers of the argument. | 

2) Integration by the method of successive approximations. f 

3) Expansion in a series in powers of small parameters entering S 

into the equation. _ | ; 

4) Expansion in a series in powers of the initial values of the t 

unknown function and its derivatives. if 

5) Method of successive approximations applied to equations for 

* 

vibrational motion. f 

6) Methods of Euler, E.L. Bravin, and others. | 

7) Method of numerical integration. 

, The first six methods do not require the use of finite differences; 

all variants of the method of numerical integration are based on the 
use of such differences. 

The principal variants of the method of numerical integration 1 

are discussed in the book by Academician A.N. Krylov, "Lectures on 
Approximate Computations", /”4_7 t 

Fundamental information on the theory of finite differences and 
on the technical features of their application to the engineering of 
artillery may be found in the book by Professor G.V. Oppokov, "Numerical 
Analysis Applied to the Science of Artillery". /~5 J g 

( Internal ballistics is an applied science possessing a 

perfectly definite technical content (the study of the motion of 
a projectile in the bore of a gun and of the laws of burning of powder) 
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and a perfectly well-defined technical objective: the cx*eation of 
Beans for plotting the curve of the speed of the projectile in"the 
bore.and for plotting the curve of the pressure of the powder 
gases in“the bore as functions of the path of the projectile and 


. _ These curves can be plotted after obtaining suitable tabulations 
giving: the functional relation between the various variable 
quantities participating in the phenomenon of a shot. The necessary 
tabulations are obtained by analyzing primary experimental data and 
those formulas which, in their simplest fora, express the relation 
existing between the initial variable quantities. 

The use of numerical analysis constitutes the subject treated 
in this section of the book.. The essence of numerical analysis, 
its specific features and its principal operations will also be 
considered here in the proper degree. _. 

CHAPTER! - NUMERICAL INTEGRATION BY FINITE DIFFERENCES 


(Written .by Professor G.V. Oppokov) 

1. APPLICATION OF NUMERICAL INTEGRATION TO THE DETERMINATION 
OF FUNCTIONS ------ 

1) Concept of Tabular Function s - t 

That variable function whose, numerical values can be chosen 
arbitrarily is usually designated as the argument or as the 
independent variable quantity, the. remaining variable quantities 
taking part~in the process under consideration are then designated 

as functions7 ...” ~ 

Let us assume that a certain independent variable quantity takes 
on a series of particular values: 
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which are separated from one another by invariable equal intervals,, 
so that: 



This interval is called the step of the argument and is 
designated by h. 

In addition to the step h, the limits of the region x Q and x R 
must be stated in the form of finite numerical values. 

Cases in which the step h is variable, or in which the variable 
x assumes infinitely large values in the region under consideration 
from Xq to x n , including the limits of the region themselves, are 
not considered at all. 

Furthermore, let some other variable quantity y also assume 
a series of particular values: 

V y r V***' *!'•••' V 


each of these particular values corresponding to one of the particular 
values of the argument x, so that it is always true_that: 


r 


y i ” ‘ 


where: 


i - 0, 1, 2 ,..n - 1, n. 
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It follows that- the functional relation between the variables 
y; and x is established in the form of a tabulation. 

, , A n example of such a relation is presented in Table 4. 

Table 4 - Pressure. Curve as a Function of Time t. 


it-:* 3 se< 


— ? — 

22 

r 

23 

23.5 

— 

24 

24.5 

25 

128 

175 

223 

274 

330 

394 

29.5 

0 

30.5 

30.75 

31 

- 

1516 

1743 

1983 

2097 

2175 



This table is the result of analysis of experimental data 
S obtained by burning in a bomb a weighed sample of strip powder grade 
JSP (1 x ,18 x 40 u) at A- 0.201. The values of t were chosen when 
evaluating, the experimental data, so that t is the argument; the 
values of pressure p were taken from the experimental curve for the 
chosen values of t. Consequently, the pressure of the powder gases 
, is a function of the time t. 

-Use has been made in Table 4 of the so-called Vertical notation, 

Which we shall adopt henceforth. In this notation the particular 
values of each of the variable quantities are invariably placed in 
one row,' while, each column contains one particular value of each of 
these variables. This notation is found to be most convenient for 
the various operations to which the functions stated by the 
1 tabulation are subjected. 

2) Finite Differences of Various Orders. 

Once a function has be^n given in table form, it is not difficult 
- "to find the so-called finite differences, of this function or, more 
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simply, its table of differences, these differences being of various 
orders. Thus, for a portion of Table 4, we obtain the following 
table of differences. 


Table 4-a. 


t-10 3 

23 

23.5 

'24 

24.5 

25 

25.5 

p 

175 

223 - 

274 

330 

3 - i 

466 

Ap 

48 

51 

56 

64 

72‘ 

- 

A 2 p 

3 

5 

8 

8 

- 

- 

A 2 p 

2 

3 

0 

- 


- 


Thus on the basis of the table of particular function values 
it is possible to compute the following differences: 

Ap 0 " P 1 " p 0 “ 223 " 175 - 48 

A Pl » P 2 - Pl - 274 - 223 - 51 

Ap 2 " P 3 " P 2 " 330 “ 274 " 56 












These new. differences are now called differences of the second 
order, or second differences. 

■ From the second differences it is also possible to compute in 
a similar.manner the third differences: 

A 3 P q - - A 2 p 0 - 5 - 3 - 2 

A 3 p x - A 2 p 2 - A 2 p 1 -8-5-3 
A 3 p 2 “ A 2 Pjj - t? p 2 -8-8-0 


followed by the fourth differences, etc. Generally, by definition, 
the k th difference equals: 

- &k ' 1 y i+ i - 

Rule . In formulating a table of differences, the number at 
the leftmost be subtracted from the number at the right in the 
same row, and the result recorded in the next lower row under the 
number at the left. 

It is useful to point out that in following this rule the 
differences with the same subscript (the function itself also being 
considered as a zero order difference) are automatically recorded 
in one and the same column. 
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Table 4-b. 


X 

x o 

X 1 

X 2 

X 3 

X 4 

y 

y o 

yi 

y 2 

y 3 

y * 

Ay 

Ay 0 

Ay i 

Ay 2 

Ay 3. 

- 

A 2 y 



A 2 y 2 


! _ 

A 3 y 

O 

>» 

% 

A 3 7! 


~ 

. 

- 


In the table of differences, in the column for i - 0, for 
example, are found differences such as: 

p 0 - 175, A Pq - 48, . ( 

a 2 p 0 - 3 - a 3 p o ' 2 

etc. 

3) Certain Properties of Finite Differences. 

1. All differences of any order can be expressed only in terms 
of particular functions of the .function itself. 

By definition of the first differences: 

Ay 0 " y l ' y 0 > iy l ’ y 2 • y l* Ay 2 • y 3 - y 2 * 

Then, by definition of the second differences: 

A \ " Ay l ' Ay 0 ' (y 2 - y i ) - <y i • V ■ H - 2y l + y o : 

A 2 yj -Ay 2 -&yj - (y 3 - y 2 >- <y 2 - - y 3 - 2 y 2 + yj 


etc. Generally: 
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&2 n ■ n+2 


- 3 >'m + >v " 

' .S ' 

The same procedure can also be applied to differences of higher 
• - ' 

orders. _- 

2. A constant number can be taken outside the difference symbol 
of any order and, conversely, can be brought inside this difference 
symbol. / 

Let: 

y - c • f(x). 

where c is a constant number. On the basis of the definition of 
differences, -we have: 

Ay - c*f (x + h) - c* f (x) - c-/~f(x + h) f (x)_7 - c* A f (x) . 

Thus, , 

A/~c*f(x)_7 - c-Af(x). 

By rewriting this relation from right to left, we obtain: 
cAf(x) -A^fcf(x)_7. 

This is -a mathematical formulation of the second part of the 
assertion under consideration. 

‘3. For an entire function of the k** 1 degree: 
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(because and h are constant numbers), which is what we set 
out to prove. 

~ 4) Determination of Coefficients of an Entire Function. 

Let us undertake the task of expressing the coefficients of. 
an entire function in terms of differences of this function. To 


do this for the case of k - 2, we shall make use of the following 


- A 2 ; Ay 4 - A 0 h 2 + Ajh; A 2 y t - 2A 0 h 2 . 


A 

0 2h2’ 


A yj . 1 t>\ + A,h, 


A 1 ' n <2A >'i - A »i>' 


Thus, finally: 



F-TS-7327-RE 


2 ^ " 

A - -i--i; 

1 2h2 


CHARLES A. MEYER& CO. Inc. 






The coefficients .for an entire third degree function can be found 
in a similar manner: 


A 


6Ay.. - 3A“y A + 2 a y i 


A 3 - y r 


These coefficients will be needed subsequently for deriving the 
formula of the interpolation function. It is of interest to note the 
fact that the coefficients of the entire function are here expressed 
in terms of differences with the same subscript i. 

The same method may be used to determine these coefficients 
in terms of differences with different subscripts. Thus, for 
example, by taking tj, t i _ 1 , t i2 > t i3 (instead of t^ t i+1 , t i+2 » 
t i+3 ), we shall find the following relations: 

* 3 yj-3 , . A 2 yi 2 + A 3 yj - 3 

0 6h 3 1 2h 2 


6 Ay. 


These relations, which the reader himself will be able to derive 
directly by using the same method, will find an application 
in the derivation of working formulas for the numerical integration 
of equations. 

5) The Pra c tical Value of Differences. 

Differences are employed for the following purposes: 
a) Determination of intermediate values of a function. 
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b) Computation of intermediate values of the argument. 

c) Factual determination of derivatives of various orders of 
a function. 

d) Determination of definite integrals. 

e.) Numerical integration of ordinary differential, equations. 

The universal character of thejte differences us'ed-in the 
operations enumerated above is reflected j.a the fact that it is 
perfectly immaterial whether the function has been stated in the 
fora of a tabulation, a diagram, or a formula. 

It has already been established that the k**V difference of an 
entire function of the k th degree is constant. It is not difficult 
to show that, conversely, if the differences of the k th order of 
a certain function are constant, the latter is an entire function 
of the k th degree. 

This proposition constitutes the cornerstone of the utilization 
of differences in all of the operations indicated above, for the 
following reason. If the formulation of the table of differences of 
a certain* function shows -that the k*k~differences are almost constant, 
we have the right to replace* our function by an entire function of 
the degree and then subject the latter function to all the 
necessary operations. ‘ 

It is this entire function which is designated as the interpolation 
■function. A general expression for it must be derived. When k - 3, 
we have: - . ' 

j - A Q (x - * a ) 3 + 'Aj/<x - XjJL + A g U - x A > + A 3 . 
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Let us introduce such a variable | that: 


x - x. + E h. 


where h is the step. Then: 


5 - 


h 


The quantity t is called the coefficient of interpolation.. 

Using the coefficient of interpolation and the general relations 
for the coefficients A^, A^, Ajj of the entire function as derived 
in Subsection 4, we shall rewrite the formula for the interpolation 
function as follows: 

e 2 6 s n* 

or, finally, after regrouping the terms: 

y - y i - -~$U - 5>A 2 y. + Mil - S)(2 - S)A 3 y r <105) • 


The interpolation function is -thus-found to be expressed in 
terms of differences of the given function up to and including those 
of the third order, with all these differences having one and the 
same subscript. 

In replacing any function (regardless of the manner in which it 
is stated) by this interpolation function,it'is categorically . 
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imperative to dir.oct attention, on the basis of the table of 
differences of the given function, to the character of the variation 
of these differences; it is necessary; and triply necessary, that 
the k** 1 differences be nearly constant. 

Obviously this condition will bo satisfied the better the smaller 
the absolute values of the differences of any order. 

" • Rule . In order to reduce the differences of any order, the• 

/ 

step must be- reduced. 

/ 

! This rule results from the following relation: 

- yOO • h^ + h* '• t- h^(y^) +6), 


where y^ is the k in -order derivative of x, and t is an infinitely 
small quantity of the highest order. The relation itself, which 
is derived at the proper stage of the differential computation, 
represents a generalization of the better-known relation between 
the increment of a function and its differential: 


vth_, 


Ay - dy + hi 

or: 

Ay - y’h + ht - h(y’ + t). L 

Thus, as the step is halved, the first differences decrease 
approximately by a factor of two, the second differences decrease 
by a factor of four, and the third differences decrease by a factor 
of eight. 
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6) Determination of Intermediate Values of a Given P uivct ion. 

To determine the intermediate values of a given function, 
regardless of how it is stated, it is necessary to formulate a 
table of differences of this function and, having made certain 
that t.he second differences are constant or are very small in 
comparison with the particular values of the function itself, replace 
it by an entire function of the second degree, i.e., by an interpolation 
function (105): 

y “ + l • - -i- |(1 - 5 )A 2 y i . (106) 

where designates the -tabular value of the given function 
corresponding to the next smaller value of the argument nearest 
to that value of x for which the intermediate value of the function 
is sought. 

Knowing the intermediate value of x for which the intermediate 
value of the given function is' sought, and having obtained from the 
table of differences of this function the nearest value of we 
can determine the value of the coefficient of interpolation £ and 
the unknown value of y by means of formula (106). 

Example 1. Given the following table of pressure impulses: 
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where the impulses are expressed in kg*dm” *sec, determine the 
pressure impulse for t - 0.020556. 

Solution. Let us formulate the table of differences: 

- • _ Table 4-d. 


t • 10 3 

20 

21 

22 

I 

Al 

A2l 

71.0 

9.6 

1.9' 

80.6 

11.5 

92.1 


Let us find the coefficient of interpolation: 


„ t - t 4 0.020556 - 0.020 

5- -- - - - 0.556. 

b 0.001 


The desired intermediate value of the function is: 

- I A +!> AI i - -U(l - ^)A 2 I i - 71.0 + 0.556 * 9.6 - —0.556(1 - 


- 0.556)1.9, 


S - 0.556; I A - 71.0; Al - 9.6 and.A 2 ^ - 1.9 


We finally obtain: 

I - 71.0 + 5.33 - 0.25 » 76 kg*dm" 2 -sec 
7) Computation of Intermediate Values of the Argument. 

Cases are sometimes encountered in practice of a contrary nature 
in which an intermediate value y of the function is given and it is 
required to find the coefficient of interpolation £ and the intermediate 
value of the argument x. Such an operation is sometimes called an' 
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inverse interpolation. . * . 

Then, discarding the last term in formula (106), we obtain as 
a first approximation: 


, s - — 
y V 

As a second approximation, we find for ? from the same formula 


Ay i “ T (1 ' ^^i 


Let us assume in the right-hand side of the above expression 
then in the second approximation: 


Ay i ~ ii (1 - 5 i )A2y i 


where 5^ is already known from the first approximation. 

Sometimes the following expression is used as the third 
approximation: 


Ay i ’ t (i - vA + t (i - v <2 - vA 
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This formula can bo easily obtained from the general expression 
for the interpolation function (105) by* placing S'outside the 
parentheses, determining a'nd substituting in the' right-hand side 
| ■ known from the second approximation. ^ 

. The desired value of x will be: 


x - + ?gh or x. - x^, + $^h. 


Example 2. We have the following table for the relative 
portion ^ of the charge: 


.It* 10 | 20 1 21 i 22 


Y 0.034 0.042 0.054 


It is desired to find the value’ of t corresponding to the inflow 
of gases 0.038. 

Solution. Let us formulate the following table of differences: 


t ■ 10 | 20 21 22 

’* H* 0.034 0.042 0.054 

0.008 0.012 

A 2 h» 0.004 - 


We obtai n as a f irst approximation: 

% 0.038 - 0.034 ! 


We find as a second approximation: 
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0.038 - 0.034 
1-T 


1 9 0.008 - \• —0.004 7 

- -(1 - t^A 2 ^ 2 


- — - 0.556. 


The desired value of the argument is: 


t - t A + ? 2 h - 0.020 + 0.556(0.024 - 0.023) - 0.020556. 


8 ) Numerical Differentiation of Functions. 

For this operation it is necessary, first of all, to formulate 
a table of differences of the given function and to make absolutely 
certain that the differences of the third order are nearly constant. 
This fact makes it possible to replace the given function (regardless 
of the manner in which it is stated) by the interpolation function 
(105): . 

y - y A + 5 • - 5 )A 2 y i + — Sd - 5 ) (2 - l)A 3 y i . 


From this we find the desired derivative: 



keeping in mind that the differences y^, Ajr^, A 2 y^, and A^y^ are 
certain constant numbers. 
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57 ’ l[ Ay i + T (25 ' “A + 7 (2 ' 6? + 3l2,A>y i] • 


E- X ' X i ; 

h 


d!_1_ 

dx h ' 


Therefore, finally: 




Such is the general working formula for numerical differentiation 
by means of differences of the giien function. It permits finding 
the value of the derivative at any value of the independent variable x 
(i.e., at any 5). This formula assumes its simplest form at 5 - 0: 


if 1 2 \ 1 

' T |* y i - 7* y M + 1* y i - 


This relation permits determining the derivative only for those 


values of the argument which appear in the table of differences of 
the given function (i.e., for 5- 0), \by applying the subscript i 
successively to each column of this taple. 
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If the stop 1) is not a prime number (for example, -), the 

738 

constant h may be introduced under, the sign of the difference of 
any desired order. Therefore, upon introducing into the process 
the auxiliary function: 


1 2 t 2 3 i o 

A$ - —Ay. ; A $ . » —A y . ; A $ » -—Ay . , 
i h '1’ i h 1 i h 


and the working formula for numerical differentiation acquires its 
final form: 


in. ] . 1 aV + i-A 3 * . 

d> 1 2 1 31 


Example 3. We have the following table for^, the relative portion 
of the charge: 

_Table 4-g._ 


t.10 3 

21 

21.5 

22.0 

22.5 

‘io 3 

42 

47 

54 

64 


t - 0.0210. 

Solution. Let us formulate the following table of differences: 
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where the limits of integration a and b are finite, this operation 
is equivalent to computing the shaded area shown in fig. 152, which 
is limited at the top by a curve representing the function whose 
integral is to be found. To compute this area Y, the usual method 
is applied first: the interval of integration (b - a) must be 
divided into n equal portions, lines normal to the O-x axis are 
then erected at the po.ints of division, and the unknown area is 
divided into elementary areas: 


The problem is then reduced to the determination of these 
elementary areas. For this purpose it is sufficient to replace the 
portion of the curve'-'f (x A > corresponding to the elementary area 
F-TS-7327-RE .635 
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by a cubic parabola, kce, 'g in mind that: 


This replacement is accomplished by means of the interpolation 


y - y* +X -Ay, - -~* (1 - ?>A 2 y,-+ Mu - ?>(2 - ?)A 3 y ., 

. i & 1 6 l 


x - x 1 + 5h, 


dx - hd?; ^ - 0; ? i+1 - 1. 


Consequently, after replacing the variables, we have: 

1 

Y i " h j - -1 5(1 - 5)A 2 yi + i-?(l - 5)(2 - |)A 3 y. d5. 

o L. J 


We integrate the right-hand side of this relation keeping in 
»ind the fact that y i# Ay i , A 2 y. i and A 3 y A , being particular values 
of the corresponding differences, .are constant numbers: 

4 - >.i v - U'4. eW. . ,•. vV’„ 
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from which, after substituting the extrame values of we obtain: 
AY, - h ^y, .+ \ Ay. - i i 2 y. + ± . 


This is the working formula employed in finding a definite 
integral by the method of numerical integration. 

To determine the unknown integral, it remains necessary 
to summate gradually and successively the individual elementary 
areas: 

Y x - Y 0 +AY q ; Y i - Y. j 4-AY^; 


Example 4. We have the following table for a function to be 
integrated; # j 


t-10 3 

8 

12 

16 

20 

p, kg*era -2 

26 

32 

48 

88 


It is desired to find the elementary area AI^, the subscript i being 
applied to the first column. \ 


Solution. Let us formulate tne following table of differences: 
_ Table W-J. _ 


t-10 3 

\ 

8 

12 

16 

P 

Ap 

6 \ 

32 

16 

40 
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10) Determination of Pressure Impulse from a Pressure- 
Bomb Test. 

As is already known from the preceding course in internal 

ballistics, the pressure impulse is expressed by the following relation: 
t 



It is clear that in order to determine it from 'ft"bomb"'test 
( yielding directly the (p, t) curve, it is necessary to compute a 

definite integral, for which the pressure of the powder gases p is 
the function being integrated and the time t is the independent 
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variable. Tins can be accomplished of course by means of one of 
the usual methods of quadratures. However, it is usually 
necessary to have a Curve for the pressure impulse as a function 
of the gas inflow^; in other words, in addition to the pressure 
impulse 1^ at the end of the burning of the powder, it also 
becomes necessary to find a series cf intermediate values for 
this quantity, which will correspond to intermediate values of the 
gas inflow. This problem is solved most simply by the method of 
numerical integration. 

Exam p le 5. Given the pressure curve (p, t) represented by 
Table 4, find the correlation (I, y) by means of the table correlating ty 
with t. ‘ 

_Table 4-k - Experimental Table of ^ as a Function of t. 


t-10 3 

0 

4 

8 

12 16 

18 20 

21 22 23 23.5 

24 24.5 

H»-10 3 

0 

2 

3 

6 14 

22 34 

42 54 79 101 

127. 155 

t-10 3 

25 

25.5 

26.0 

26.5 27.0 

27.5 28.0 

28.5 29.0 29.5 30.0 

30.5 31.0 

4>-10 3 

186 

221- 

260 

304 354 

409 470 

539 619 711 813 

919 1000 


Solution. We perform all computations on the working form in 
Table 5, without considering third differences, and designating: 

1 12 1 i o 

£« y + — Ay - —A y - p + —Ap - — A p: 

2 * 12 2 12 

from which it follows that: 


i - AY i - h k - “ ^ 2y i " h * r ' 


F-TS-7327-RE 


■CHARLES A. MCYCR& CO. Inc. 








175 223 



Ul « hr 18.0 21.4 25.2 29.5 - 34.4 39.8 46.1 53.1 61.2 70.6 81.4 93.4 - I- fe 1 































Table 5 - Computation of Pressure Impula 









The computations in Table 5 are conducted in rows, proceeding 
from left to right in each row and downward from one row to the 
next, in the following manner. 

3 

1) The first row is filled with values of the argument, t*10 . 

2) The numbers for the pressure are taken from Table 4 (page 
616). - 

3) From each number in the second row there is subtracted 
the preceding number in the same row, and the result is written 
under the left-hand number of this pair: 

23 - 21 - 2; 26 - 23 - 3; etc. 

4) The numbers in the' fourth row are obtained in the same 
manner as in the preceding row: 

3-2-1; 6-3- 3; etc. 

5) The numbers^of the second row are repeated. 

6) The numbers in the third row are halved. 

7) The numbers in the fourth row are divided by 12 and written 
with the opposite sign. 

8) - The^numbers in the three, preceding rows are added. 

9) The \numbers in the .eighth row are multiplied by the step h, 

the decimal point being correctly placed (cf. Example 4) to convert 

_ \ „ 

kg*sec*cm~ z into kg*sec*dm . 

10) To obtain the next succeeding value of 1^,. it is necessary 
to add to the preceding value of the corresponding increment 

M i-1 : 

“ x i-l + Al i-1» 
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or, in other words, it is necessary tc add tlie two numbers in the 
ninth and tenth ro-s of the preceding column: 

0 *■ 8.8 - 8 . 8 ; 

8.8 + 10.0 - 18.8,etc. 

By taking the values for the gas inflow and the pressure 
impulse 1 corresponding to the same instants t, it is easy to 
obtain the desired correlation between 1 and ^(cf. Table 6). 


Table 6 - Correlation between I and ^ for Te st. 


I kg*sec*dm -2 

0 

& 

19 

30 

45 

56 

71 

81 

92 

107 

117 

129 

144 162 

vy-10 3 . 

0 

2 

3 

6 

14 

22 

34 

42 

54 

79 

101 

127 

155 186 

I kg*sec'dm" 2 

184 

209 

238 

273 

313 

359 

412 

473 

544 

625 

718 

770 

823 

CO 

o 

il 

221 

260 

304 

354 

409 

470 

539 

613 

711 

813 

919 

967 

1000 


Table 6 will be needed subsequently for the solution of the 
principal problem of internal ballistics (determination of the curves 
for the speed of the projectile and for the pressure of the powder 
gases as a function of the path of the projectile). This solution 
must be prepared by a discussion of the necessary theory,' which will 
be undertaken in the next section. 

• 2. NUMERICAL INTEGRATION OF DIFFERENTIAL EQUATIONS 
11) Numerical Integration of First-Order Equation 
We have such an equation in its general form: 

F(y, y’, x) - 0. 
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It is composed of the following elements: the independent 
vari:>ble x, the unknown function y, and its fjrst derivative y'. 


•It is necessary, first of all, to determine this derivative from 
the general equation: 


y' - f(y, x). 


It is this equation which is solved by the method of numerical 
integration, regardless of the form of the function which constitutes 
its right-hand side. In practice, what is usually found is not 
the total integral of the last equation (in which case, of course, 
numerical integration cannot be employed), but, rather, a partial 
integral, in consequence of which it is necessary to indicate the 
’’initial conditions", i.e., the values of y^ and x q . In addition, 
the value of x r at the end of the finite region of integration must 
also be known. Under these conditions the desired function will be: 



x 0 


or 


y 


f(y, x)dx. 


The problem is consequently reduced to the finding of a definite 
integral whose characteristic feature is such that the part of the 
function to be integrated is played by the derivative of the function 
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•It would seem that, as on the previous occasion involving the finding 
of a definite integral, it might be possible to perform computations 
by rows and to make use of a working formula containing differences 
of an auxiliary function with they4;.mc subscript i. In actual 
practice,, however, this procedure cannot be adopted. As a matter 
of fact, the last formula shows that the function to be integrated 
includes the original function itself, and we ttyus obtain a "vicious 
circle;" in order to determine the particular value' y x of the function, 
we must know the particular value 4>^, and in c/rder to determine 
this particular value of the auxiliary function: 

*i “ f ^i> x i )h 

we must have the value of y^ which is a constituent of- the auxiliary 
function. 

It is now clear that in performing numerical integration of the 
equations it is not possible to perform computations by rows, it" 
being necessary instead to prcceed gradually, step by step, performing 
the required computations downward along each column and to the 
right from one column to the next. 

Let us further assume that the values and ^ 

of the auxiliary function have already been found. With their aid 
it is possible to compute the following differences. 
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$ 

0 • 

4> . _ 

4> 



i-3 

i-2 

l-l 


A4> 

A *i-3 

A4, i-2 



aH 

*Vs 

* 2 *i-2 



a\ 






mi i + 4 a *j - 


in the table. Consequently, it becomes necessary to derive an 
additional formula containing the differences: 

; A \. 2 ’ A-3' 

already contained in'the table. 

For this purpose, as on the previous occasion, we shall replace 
the functloh to be integrated, y’ - f(y, x), by the following 
interpolation function: 

3 2 

y’’ » A (x «* 3^) + A x {x ~ x A ) + A 2 (x - x ^ + Ag, 
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The above table shows that in order to determine Ay^, it is not 
possible to make use cf the following formula for determining definite 
integrals: 


yz*-.,.,.. 









f(y, x)dx 


4y i 


i - J f(y, x)dx - J /~*Aq(x - x A ) 3 + Aj(x - Xj ) 2 + A 2 (x - x A ) + A^/dx - 


—A»(x - x .) 4 + -i-Ai(x _ x ) 3 .+ — A (x - x ,) 2 + Ao(x - x.) 
4 0 1 3 1 1 22 1 *> 1 


whence 

1,4 1 . u 3 1 . J2 

Ay i " T V + T A l h + “2 A 2 h + A 3 h> 

because 

X 1 + 1 -x. -h. 


In'contrast with the preceding case, we shall substitute here 

i 

the following values for the coefficients: 








Ayj - My! + -i-Ayi.! + ^y'i-2 + i i3y ;-3 K 


It is usually convenient to make use of the following auxiliary 
function: 


$ - hy 




the derivative y' being computed in accordance with the given equation: 


y* - f (x, y). 


Then: 

hy’ -$ 4 ; hA 2 y i -2 m * 2 *X-2 ; h * y i-3 " A ^i-3» 

and the working formula for the numerical integration of ordinary 
differential equations of the first order will assume the following 
form: 

A y 4 - + ^a\_ 2 + T fl3 *i-3'. 

This formula now includes differences of the auxiliary function 4? 
provided with precisely those subscripts which are required in 
accordance with the table of difference presented above. 

It will be seen from the same table that in order to determine 
Ay^ it becomes necessary to take differences from different columns, 
which, of course, is inconvenient. For this’reason it is desirable 
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to modify the system of writing the differences in such n manner as 
to have in the same column differences with different subscripts, 
na me1y ; 


Table 6-b. 


$> ■ 

*i-3 

'*1-2 

Vi 


A* 


A *i-3 

A *i-2 

> 

•S' 

1 

aH 



a2 +1-3 

A 2 *i _ 2 

A 3 4> 






In such a procedure each column will contain those differences 
which are necessary for the application of the new working formula 
for numerical integration. 

Rule f or writ in g new differences. From the right-hand number 
in a row there is subtracted the left-hand number in the same row, 
but the result is written under the right-hand number. 

Finally, from the last table of differences there follows the 
.most disagreeable characteristic feature of the numerical integration 
of equations: in order to determine Ay^ it is necessary to know the 
differences: 


A*. 


A 


i.e., to have the following particular values of the function: 


These will become known in the course of the process; but at 
the start of the integration only one particular value of this function, 
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namely^, is Known from the originally stated equation for y*. 
In order to find- % 



it is necessary to employ special supplementary methods, of which 
the .most commonly used is the method cf successive approximations. 

The essence of this method resides in the fact that, 9 .t the start 
of the computations, there is adopted a gradual and /tepwise advance, 
each new approximation permitting one additional step, in which 

those differences that have already appeared earlier are utilized. 

/ 

There exist several variants cf this method, one of which, the 
most exact one, can be best studied by the aid of a concrete example. 

Example 6 . Solve by the method of numerical integration the 
equation: 


y' - x - y 

over the interval 0 ; 0»5, with~a stop of h ■ 0 . 1 , if y - 1 when 

x - 0. 

Solution. All the computations are presented in Table 7. 
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Integration of Equation 































Tabla 7 - Integration of Equation J 



Second ‘ • 

" approximation Third approximation Normal computations 





























First Approximation 


First Column 

13) This is the leading row and is filled first; in the given 
case y 0 •» 1.0000 from the conditions of the problem. 

‘ 1 ) Xq - 0 from the conditions of the problem. 

2) -yg - -1.0000 in accordance with the thirteenth row. 

3) yg - -1; in accordance with the given equation for y', 
we add the numbers in the first and second rows above. 

4) - -0.1000; the number in the third row is multiplied' 
by the step h « 0.1, since: 

*0 ■ h> 0 

The fifth, sixth and seventh rows are not filled in for lack 
of data. Then: 

8 ) $q - -0 .1000; the number in the fourth row is repeated. 

The ninth, tenth and eleventh rows are not filled in for lack 
of data. Thereupon: 

12) A yg - -0.1000; in this case it is necessary to repeat 
the number in the eighth row. 

Second Column (of the First Approximation ) 

13) Leading row; y^ - 0.9000; we add (algebraically) tbe numbers 
in the twelfth and thirteenth rows of the first column, since: 

y-1 - y 0 + Ay 0 ; 

1) Xj - 0.1; we add the step h. 
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.2) -y^ - -0.9000 in accordance with the thirteenth row. 

3) y^ « -0.8000; in accordance with the given equation for y' 

we add the numbers in the first and second rows of this column (above). 

4) - -0.0800; the number in the third row is multiplied 
by the step h - 0.1, 

Thus a number appears in the fifth row; - 200, for which 

we subtract from the number (-0.0800) in the fourth row above the 
number (-1.0000) at the left in the same row, omitting the zeros 
at the left for greater convenience. 

Vedge of First Approximation. 

It consists merely of a single number and is surrounded by 

a heavy line. . To obtain it, we postulate that 

- 200 . 

Second Approximation 

First Column . 

13) Leading row: - 1; 1) - 0. 

The second and third rows are not filled, as this is no longer 
necessary (cf. first column of first approximation). Ihen: 

4) $ 0 - -0. 1000, as in the same row of the first approximation. 

- Omitting the fifth, sixth, and seventh rows, we write: 

8) <t> 0 - -0.1000; the number in the fourth row is repeated. 

9) — • 100; we take the difference from the wedge 

2-1 -1 
of the first approximation. 

The tenth and eleventh rows are not filled for lack of necessary 

data. 
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12) Ay 0 - -O.ODOO; we add the numbers in the eighth and ninth 
rows above. 

The second column of this approximation is formed in an 
analogous manner, with the only difference that numbers appear 
in the first and second rows (cf. the number in the thirteenth row 
of this column) and in the fifth row, for which we subtract from 
the number (-0.0810) in th£ fourth row the number (-0.1000) in the 
same row, but on the left,- The number in the leading row is obtained 
by adding the numbers in /the twelfth and thirteenth rows of the 
preceding column. 

The third column is distinguished from second in that a number 
(-17) appears in the sixth row, for which we subtract from the number 
173 in the fifth row the number 190 at the left, and another number 

(-8) appears in the tenth row, for which the number in the sixth 
5 10 

row is multiplied by — - —. 

12 24 

Filling in the Wedge of the Second Approximation . 

In the sixth line we repeat twice the number (-18), noting that: 

A 2 $ - A 2 $ • - a 2 * - -is. 

-2 -1 0 

The number 208 in the fifth row of the first column of the second 
approximation will appear in accordance with the definition of the 
second differences: 

-A*> 0 - A4>_ x , 

from which: -■“'V / 
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- A$>_ 0 - - 11*0 - (-18) - 208. 

Th ird Approx»mo tion. 

In the first column, after the thirteenth and first rows are 
filled, the places in the second, third, fifth, sixth, and seventh 
rovs are left blank. Then: 

8) The number in the fourth row is repeated. 

0) A$ j " 104; the number 208 in the wedge of the second 

approximation is halved. 

5 2 

10) —A # 0 : " -8; the number (-18) in the wedge of the second 

12 "* 

approximation is multiplied by 

The eleventh row is omitted. Thereupon: 

12) The numbers in the eighth, ninth, and tenth rows above are 

added: 


-0.1000 + 0.0104 - 0.0008 - -0.0804. 

Id the second column, after the thirteenth, first, second, 
third, fourth and fifth rows are filled, the sixth and seventh rows 
are omitted. Then: 

8 ) The Dumber in the fourth row is repeated. 

9 ) — A4 a • 85; the number in the fifth row of this column 
2 0 

is halved. 

10) _ • -8; the Dumber (-18) ia the secood column of the 

12 -2 

preceding wedge is multiplied by 

The eleventh row is likewise omitted. Thereupon: 
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12) The numbers in the eighth, ninth, and tenth rows above are 
added. 

In the third column, only the thirteenth, first, second, third, 
fourth, fifth and sixth rows are filled. There appears for the 
first time a number in the seventh row: 

-18 - (-17) - -1. 

Filling in the Wed g e of the Th i rd Approximation. 

We make: 

* 3 *-3 -^*-2 “AVI "A, - - 1 * 

The number (-16) in the sixth row of the second column of the 
wedge will appear from the definition of the third differences: 

A 2 $_l - A 3 *..! - “17 -(-1) - -16. 

In an analogous manner, the blanks in the rows of the first 
column of the wedge will be filled as follows: 

A 2 *_ 2 “A 2 *,! ~aV 2 “ ~ 16 " ( “ 1) * 

•> 2 

- A* 0 - A - lyO- (-16) - 206. 

For convenience of computation, the numbers of this wedge are 
transferred to their places in the columns of the subsequent normal 
computations. The approximations are considered to be completed 
because of the cioseness of the results of the third and second 
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approximations, respectively. 

Normal Computations. 

First Co lumn. 

13) As before, y 0 - 1.0000; 2) x^ - 0, 

The second and third rows need not be filled. The numbers in the 
fifth, sixth, and seventh rows are already in place. Thereupon: 

8) - -0. 1000; the number in the fourth row is repeated. 

9) —A* , • 103; the number 203 in the fifth row is halved. 

2 -1 

10) —a 2 * o • -6; the number (-15) in the sixth row is 

12 in 
multiplied by — . 

11) —A 2 * , - 0; the number (-1) in the seventh^rbw is 
multiplied by 

12) Ay^ - -0.0003; the numbers in the eighth, ninth, tenth, and 
eleventh rows are added. 

The subsequent columns are filled in exactly .he same manner. 

Attention is once again directed to the crder of computation in 
each column of normal computations. 

a) First of all, the leading thirteenth row is filled by adding 
the numbers in the twelfth and thirteenth rows of the preceding 
column, because 

* *i-i +A *i-i ; 


b) Thereupon, the spaces in each column are filled from top 
to bottom without omissions. 
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12) Use ot Numerical I 'itepra I ion oj the F 1 r st-Order Fq u.i t ion 
frith Argumejnt V. 

The method of numerical integration must be applied to the 
solution of,the following equation: 


fo^l 


dl 

dv 


1 

v 2 

n P 


-<L + />, 


f- 




In the physical law of pocder burning the correlation between 
and v is given by a table, whereas in the geometric law: 


% 


K6 0 «?m 

i- —-— v + --—v 

s. k 


Since it is also necessary as a rule to find the pressure curve, 
it is preferable to make use of the following equation: 

•p-v £ - p B> 

whence 

!L • t? 2 . 

dv s p * 
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•r.frT^.-.aaans; 




where 


2 


fu» _np 


It is not difficult to see that the auxiliary function ♦ is 
found in the given case from the following relations: 


'♦‘‘a - <a<K! P- - 


8 l y + ( 


As for the working forsula of nuserical integration, it has the usual 


, » 1 5 2, 3 3 

A, i + + n A *1-2 + t a *i- 3 - 


The purpose of the preliminary computations is to determine all the 


t* - 32 ; 3) „» . 2 -L 

* «« "p « 


« 7-T ‘ 5 > ‘.’T “-T : 


f 1 

~ + a - — 
P n S 
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r. 


I 


To determine the step oi integration h, the speed of the projectile| 
vj,- at the end of powder burning must be known, and then: 


where n - the number of sections - is taken in the range of 10-40, 
depending upon the density of loading (as this density increases, 
the number n must also be increased). 

While determining | and p, the same working form may be used 
to find the time of motion of the projectile in accordance with the 
following relation: 


dt 


. *m 


dv 


in the form of a definite integral. 

Example 7. Determine the projectile velocity and gas pressure 
curves for a 76 mm gun, given the following conditions: - 1.654; 

s - 0.4693; / A - 18.44; q - 6.5; p Q - 30,000; f - 900,000; a - 1; 

8 - 1.6; 0 - 0.2; 0.930; «f • 1.05; g - 98.1, using strip powder 


type SP (1 x 18 x 40 mm). It is assumed that the powder burns in 
conformity with the physical law derived from the pressure-bomb test 
presented in Table 5 (correlation between I and^). 

To start with, we find: 

log ~ - 0.2970; log ^ - T.1709; log vj - log 8.0803; 


v 2 - log ~ \ " 8 -° 

n P M «f ® q / 


1 o ( 1 - y) • 2 * 284 ; lo « “ 1.8710; - 0.038; log ^ - 6.2512. 
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Fro* the value found we compute the value of t y by inverse 
Interpolation with the aid of the data in Table 2 (cf. Example 2, 
Subsection 7) 

tg - 0.020556 sec 

For this value of t Q we determine the initial impulse 1^ by 
direct Interpolation with the aid of Table 4 (cf. Example 1, 
Subsection 6): 

I 0 - 76 kg*sec*dm” 2 

We find the velocity of the projectile at the end of burning 
of the powder, keeping in mind that - 823 (cf. Table 5): 

V - — (I - I A ) - 5040 dm.sec" 1 
K ** K 0 

of sections, we obtain the 


- 6.5723. 

the ty, I) curve to the 
for I, 1 mm - 1 kg*sec•dm” 2 . 


Bt choosing n - 20 as the number 
step as being: 

V K -l 

h - - 252 dm* sec 1 

20 


so that 



On the basis of Table 5, we plot 
following scale: for 4*, 1 mm - 0.001; 


CM-* 'I* 1- *• V-VuK 4 vO 
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We road on it tin* values of the gas inflow corresponding to t> ,ual 
intervals 

s 

8211 - 76 , , _o 

- - J7.J5 kg*see*dm “ 

20 

foi the pressure impulse, or, what is the same thing, to e^ual 
intervals h - 252 dm*sec -1 for the velocity of the projectile. 

These values of v^'are given subsequently in the working form 
used for performing the computations (in the fifth row). 
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6.2008 


-0.002 

-0.005 

0.164 

0.222 

1.8710 

1.8710 

1.2201 

1.3560 

1.0911 

1.2270 


0.262 

0.385_ 

0.573_ 

2.423 

2.500 

J2.646 

6.2512 

6.2512 

6.2512 

1.2148 

1.3464 

1.4393 

1.6257 

1.6021 

1.5774 

5.0817 

5.1997 

5.2679 


1.5723 

1.5723 

. 1.5723 

1.5723 

2.7024 

2.8785 

1.0034 

3.1004 

6.9183 

_ 6.8003 _ 

6.7321 

_6^6850 

1.1930 

1.2511 

1.3078 

1.3577 

0.156 

0.178 

1 ^ 0.203 

0.228 

24 

22 

25 

25 












































































































the working form can be broken down into three main sections. 

The upper section contains lines 1-24 and is reserved for computations 
necessary for determining the logarithm of the auxiliary function 
in accordance with the relations (*). Consequently, the pressure 
Of the powder gases is also found at the same time. The values for 
the gas inflow in line 5 are first read from the ty, I) curve in 
the case of the physical law of burning of-the powder or computed 
in accordance with the corresponding formula in the case of the 
geometric law of burning. 

The upper section is the most involved part of the work and 
aakes the use of four-place logarithms obligatory. 

In the middle section, which consists of lines 25-34, are 
performed computations necessary for the use of the working formula 
for the numerical Integration of a first-order equation: 


— A 2 *. + — aV 

12 i-2 8 i-3 v 


and for finding the path of the projectile: ^ +A *i-1‘ 

Finally, the lower section, consisting of lines 35-47, is reserved 
for finding the time of motion of the projectile in the bore in the 
form of a definite integral, using the following working formula: 


At, - F, ♦ -i-AF, - ± A\ ♦ ±L\. 


\ 

. This section is best filled after completion of the numerical 
integration of the equation for the path of the projectile: 
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i.e., after the two sections above it are filled. 

It is useful to point out some of the characteristic features 
of the computations. 

' a) The path of the projectile is determined with an accuracy 
of 0.001 dm; the time t, within 0.00001 seconds. 

b) To facilitate computations, only one approximation is 
made, and only one wedge is filled out (cf. line 26, first column 
of the working form). 

c) At first, the differences of the auxiliary function up to 
and including the second order are introduced. 

d) The third differences of this function are utilized only 
after p^ is passed. 

e) The second difference of the function $ appears only in 
the second column of the normal computations. 

f) The rows which are not filled are: lines 2, 3, 4, 6, 7, 21, 

22, 23, 24, 27, 28, 30, 31, 32, and 35-47 in the first column, lines 
27-33 and 35-47 in the second column of the approximation, lines 2-24, 

26, 27, 28, 31, and 32 of the first column of normal computations, 
and, finally, lines 2-12, 27, 28, 31, and 32 of the second column of 
normal computations. 

g) Line 34 is the leading row and is filled first in each column. 

For a more definite conception of the character of the work Involved, 

when filling in the two upper sections of the working form it would 
be desirable to describe in detail the computations performed for 
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OHO of trie columns, for example, for the column of normal computations 
-iut * - 504. 

34) The ,;•«?»burs ’n lines 33 and 34 of the preceding column are 

added: 


•1 


l - 0.198 + 


0.064 - 0.26?. 


^i~I + 


*( 


1) Ve add the integration step h: 

v - 252 + h - 252 + 252 - 504; 

2) The number in the preceding row, v - 504, gives: 

2 log v - 5.4048; 

3) The complement of the logarithm of is taken: 


-log v„ - 9.9197; 

Op 

4) The logarithms in lines 2 and 3 are added: 

g jL - 5.4048 +*9.9197 -'3.3245; 

V °P 

5) The value of ^ is taken from the (^, I ) curvex 


- 0.166; 
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~.**“i*#*$m 

6) The logarithm in line 4 is used to find the number vithin 

O.Obl: . 

✓ 

v 2 

- 0 . 002 , 

V n P 

7) The number in line 6 is subtracted from the number in line 
5: 

v 2 

-- 0.166 - 0.002 - 0.164; 

v 2 

np 

8) The logarithm of has been found by preliminary computations: 

log l a - 1.8710; 

9) The number in line 5 gives the logarithm: 

log - 7.2201; 

10) Ihe logarithms in lines 8 and 9 above are added: 
log l a " T. 8710 +7.2201 - 1.0911; 

11) The reduced length l^ has been found by preliminary computations: 

L - 2.284; 

A 

12) The logarithm in line 10 is used to find the number vithin 

0 . 001 : 
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ft - 0.123; 


13) The number in line 12 ts subtracted from the number in 

line 11: 


" 2.284 - 0.123 - 2.161; 


14/ The number in the leading row, line 34, is taken: 

/ I - 0.262, 

15) The numbers in lines 13 and 14 are added: 


l^+ t - 2.161 + 0.262 - 2.423; 


16) This logarithm has been obtained by preliminary computations 


log — - 6.2512; 


17) The logarithm of the number in line 7 is found: 

“■Mb'"" 


18) The complement of the logarithm of the number in line 15 

aken: 


-lo*(l^ + 1) - -log 2.423 - 1.6257; 
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19) The'logarithms in lines 16, 17, and 18 are added: 


log p - 6.2512 + 1.2148 + 1.6257 - 5.0817; 

20) The logarithm in line 19 is used to find the pressure 
at- v • 504 : 

p - 1207 kg*cm"2 

21) This logarithm has been obtained by preliminary computations: 



22) Prom the number in line 1 we have: 

log v - 2.7024 

23) The complement of the logarithm in line 19 is taken: 

-log p - -5.0817 - 6.9183 

24) The logarithms in lines 21, 22, and 23 are added: 
log 4> - 1.5723 + 2.7024 +*6.9183 - 1.1930. 

We now proceed to the next, middle, section. \ 

25) The logarithm in line 24 is used to obtain the number with 

an accuracy of 0.001: \ 

* 2 - 0.156! ^ 
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It is now possible to take the next step, i.e. f proceed to 
the next column with v - 756, starting the computations therein, 

.as always, fro* the leading row, line 34. 

At the muzzle^we have: - 583^m/sec; p^ - 625 kg'em 2 . 

13) Numerical Integration of the Second-Order Differential 
Equation. 

In the general case, a differential equation of the second 
order contains the following components: the independent variable x, 
the function y itself, the first derivative y* of the function y 
with respect to x, and the secopd derivative y" of the same function 
with respect to x. Consequently, this equation can be represented 
in the following for*: 

*(x, y» y*# y") - o; 

where the symbol F represents an elementary function or a combination 
of such elementary functions. 

Numerical integration makes it possible to solve this equation 
in any form, provided only that the given equation permits the 
deterslnatlon of the second derivative as an explicit function of 
all the remaining variable quantities: 

y" - *(*, y# y’)* 

In this numerical integration the increment of the derivative 
Ay‘ is determined by the aid of the following formula: 
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Ay^ - h 


or 








3 3 


A *. 


where use is made of the following new auxiliary function: 
4> - hy", 


expressed in terms of the second derivative of the desired function. 
It is not necessary io derive the fundamental formula for Ay; by 
means of y" and differences of this second derivative, inasmuch 
as this formula is already obtained from the previously derived 
relation 


Ay. - h y 


i + T y [.i 


5 2 , 33 . 

+ 12 A y l-2 + T* K- 


by simply replacing y^ by y; and y; by y^’. 

On the other hand it is necessary to derive a working formula 
for the increment Ay of the function itself in terms of differences 
of the second derivative rather than of the first. Such a formula 
should simplify the work, because in computing Ay it will be possible 
to make use of the already available differences of the auxiliary 
function <t>, and it will not be necessary to find the differences 
of still another auxiliary function: 
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For this purpose, we shall utilize the equality 

x r 

Ay A - \ y'dx. 


But in all cases 

/ y' - y 4 * + Ay 


I - 


i 




y"dx. 


Let us replace here the derivative y" by the following inter¬ 
polation function: 

3 2 

y" • Aq(x - X A > + Aj(x - x i ) + A 2 (x - x A ) + A 3> . 


5,1 ' y i + 


A 

j" Z r A Q (x - x,) 3 + A a (x - x i ) 2 + A 2 (x - x A > + A 3 _7 dx - 


y i + T A 0 (x “ x i )4 + A l (x ’ x i )3 + \ A 2 (x “ x i )2 + A 3 (x “ x i ) 


that, consequently 
x i+i 

.( AT 


'I 


y i + T A o (x “ x i )4 4 T A 1 (X ” x i )3 + 
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~- A 2 (x - x i > 2 + A 3 (x - Xj) 


dx - y’(x - x A ) •. ± A 0 (x - x 4 ) 5 


X i+1 

+ T2 A 1 (X “ X i )4 + 7 A 2 <X " X i )3 + 7 A 3 (X “ x i )2 

x i 


Ay ■ h y! + i Ah 4 + 
i i 20 0 12 


i 1 2 1 

Ah 3 + — A_h + —A_h 
1 6 2 2 3 


x^i - x, - h. 


We shall substitute into the last expression the following values 


for the coefficients: 

y i’-3. 


A; 


A ° " 6h3 


Jy't-2 * * 3 yi'-3. 


6Ayf. 1 + 3A 2 yJ'_ 2 + 2A 3 yJ'_ 3 


h - yj'- 


because the part of the function is here played by the derivative 

1 

y", which has been replaced by the interpolation function. 

After obvious transformations this substitution gives: 
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atj -»y;+ i A »j-i + i'fy-z + rss A3y ;'-: 


Ay A - h \ y [ + h 


/ 1 . 1 .. 1 ,2 „ 1 a 3 „ 

\* y i + 6 Ay i-1 + T A y i-2 + IS A V 


if, instead of the iuconvenient number we take the closely 
spproching it number: 

18 _ 1 
180 10 


It remains necessary to make use once again of the auxiliary 


and then, finally, in the numerical integration of the second- 
order differential equation, we shall have to deal with working 
formulas of the following type: 

Ay i * *1 + T 4 *l-1 + H * 2 *i-2 + 4 a3 *1-3 ! 


T 1 - h ( y i + i*i + + 7*1-2 + 


To simplify the computations, it is better, however, to conduct 
the numerical integration with so small a step h as would enable us 
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to do without third differences. 



At the start of the numerical integration, as before, most 
frequent use is made of the method of successive approximations. 

The number of rows in the lower part of the form will increase 
(in comparison with the number of rows in the lower part of the form 
lor the integration of the/first-order equation) by at least seven 
rows, because in addition to the ten rows: 

♦ ,A*,A 2 *.A 3 *, rfa 2 *. J^.r, Ay' - hr, y. 


corresponding to the lower part of the form for the numerical 
integration of the first-order equation, there will also appear 
the following additional rows: 


4*’ ‘ hr i’ y - 


In the method under consideration, the computations themselves 
are in no way different from similar computations for the solution 
of the first-order equation, so that the need for citing a special 
example is obviated. 

14) Use of Numerical Integration of Second-Order 
Equations with Argument t. 

In this process the leading part is played by the equation for 
the forward motion of the projectile: 


£ - sp, 


from which it follows that 
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if the phenomenon of recoil is not considered in the explicit form. 

Equation (107) is subject to numerical integration by means 

of the following relations: 

2 

♦ V 

p-T-rrr - 


which have already been'employed earlier. The auxiliary function $ 
is equal in this case to: 

♦ - hi» - - P . 

t «pm 

where the step h is already a certain time interval of the motion 
of the projectile in the bore. This step is chosen in advance. 

The working formulas for the numerical integration will be: 

1 -5 ‘2 3 3 

i i 2 1-1 12 1-2 8 l ~ 3 

A 'i ’ h ( v i + T*i + T 4 *i-i + tA - 2 + To* 3 * 1-3)' 

because the function y to be determined represents in this case the 
relative path / of the projectile in the bore, and its derivative y’ 
Is the velocity of the projectile y.. 
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The preliminary computations are reduced to the determination 
of constants: 


1) log ~ ; 2) log 3) log ; 


4)/^ ; 5) log t Q ; H) q» Q ; 


7) log 


as in the case of numerical integration of a first-order equation. 


I^r 


The gas inflow 4> is read off the (I, curve, but, in contrast 
with the preceding case, it is necessary here to find for each point, 
during the process of integration itself: 


, *Pm 


because the integration gives values of the velocity v which are 
separated from one another by unequal intervals. Once the value of 
the pressure impulse is had, we can obtain the required value of 
from the I, curve. Consequently, this curve cannot be utilized 
in advance, but it must be available in the course of the entire work 
of numerical integration. 

In the case of the geometric law of burning, the I, ^ curve is 
replaced by the following relation: 

_2 


+ "+0 
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Fig. 153 - Velocity, Path, and Pressure Curves. 

1) Tine; 2) pressure curve; 3) velocity curve; 

4) path curve. 

Figure 153 contains curves for p, v, and l as functions of time 
t, obtained by numerical integration vith respect to the argument t. 
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Co»puta“,on Tc 




















Computa 

.on Form 

with Arw 

25 

50 

0 

1.1700 * 

1.74U4 

1.1700 

2.0755 


0.0113 

1.2464 

























































Compu t a Ion Fo r « with - Ar^uiM -nt t. 
5»5 1 5ft I! ft ... oT 1 


1.1709 1.1709 
I.7404 2.0755 
0.9113 1.2464 


1.170jI 1.1709 
2.0 8641 2.3160 
1.2575 I J. 4869 



3.4808 4.1510 
9.9197 9.9197 
5.4005 4.0707 

0.045 0.057 



1.8710 1.87i0 
2.6532 2.7559 
2.5242 2.6269 

2.284 2.284 
-v.033 -U.U42 


2.251 

0.007 

~ 2.242 

0.029 

2.258 

2.271 

6.2512 

6.2512 

2.6532 

2 .7559 

1.6463 

1.6438 

4.5509 

4.6509 

3.2370 ' 

•9.-*270 




153 

5.0422 I 5.4304, 
9.9197 9.9197 



1.8710 

1.8710 

1.8710 

2.8976 

1.0755 

1.2355 

2.7686 

2.9465 

1.1065 

2.284 

2.284 

2.284 

-o.oo» 

-U.OM« 

-U.12M 

2.225 

2.195 

2.156 

0.069 

0.135 

0.239 

2.294 

2.330 

2.395 

6.2512 

6.Ql2 

6.2512 

2.8976 

1.0719 

1.2304 

1.6394 

l.Cd26 

1.6209 

4.7882 

4.9557 

5.1025 

3.2270 

3.2270 

3.2270 

— 

1. _ 
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CHAPTER 2. SOLUTION BY EXPANS I ON IN TAYLOR’S SKKI KS 

Integration of equations of internal ballistics describing the 
relations existing between^ the fundamental elements of a shot leads 
to rather complex integrals, which can be solved only by means of y 
quadratures with any desired degree of accuracy (Professor Drozdov's 
solution). In this connection, in order to make integration possible, 
some of the variable parameters in the fundamental equations (9, u^, 

Pq, etc.) are usually assumed to be constant. 

The method of numerical integration makes it possible to solve 
a system of differential equations not only without simplifying the 
functions entering same, but even by assigning variable values to 
those quantities which are usually assumed to be constant. This 
makes It possible to solve the problem on the basis of any desired 
hypotheses concerning the character of burning of the powder, the 
law of resistance to motion, the design of the bore, etc. 

In solving the problem by the method of numerical integration, 
it is necessary to proceed successively from one value of the 
argument to another by the addition of its finite differences, starting 
at the very beginning. For this reason, for example, it is not 
possible to determine in advance the values of p ffl or the values of 
the variables v, ( and p corresponding to the end of burning of the 
powder, it being net* ary, instead, to compute successively, point 
by point, the elements of the curves of pressure p, the path of the 
projectile I, its velocity v, etc. This constitutes the disadvantage 
of this method. Moreover, the method of numerical integration gives 
the relation between the Individual variables only in the form of 
numerical tables, rather than in the form of analytical formulas. 
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In spite of these disadvantages, the method of numerical 
integration may serve .is a means for indirect verification of the 
degree of accuracy obtained by the aid of the various approximate 
analytical methods available. 

In this connection, when developing new theoretical piobletns, 
numerical integration may bo utilized for determining the errors 
involved in the transition from exact equations and relations to 
others that arc less exact but more convenient from the analytical 
point of view. Numerical integration may be employed with equal 
success both in the case of the geometric law of burning and in 
the case of the more complex physical law of burning; it may also 
be applied to barrels having a bore of variable cross section. 

In the USSR the method of numerical integration was first 
applied to the solution of ballistic problems by V.M. Trofimov in 
1918. This method has been developed in great detail by G.V. 

Oppokov, who employed the method of finite differences (1931-1938) 
discussed above. 

In 1932 P.V. Melentyev proposed to apply the method of expansion 
in laylor's series for the numerical solution of equations in 
ballistics, and this method, after being subjected to certain 
modifications, has been found to be sufficiently convenient. 

Investigation of the fundamental relations expressing the 
conditions accompanying a shot shows that all the principal elements 
(/, v, ^, and p) can be expressed in one fora or another as functions 
of the path I and of its derivatives with respect to time up to the 
third order inclusive. 


S 





CHAM;.S' 4. V'TiT.R & CO- Sr.c. 




F-TS-7327-RE 


681 






As a natter of fact, taking the fundamental system of equations 
expressing the relationship between the elements of a shot, we obtain 


the following. 

1) The fundamental equation of pyrodynamics or the equation 
of transformation of energy: 

1) - fw* - ~ «fmv 2 . 

2) The equations expressing the law of burning of the powder: 


^ - Xz(l + Xz); 


de 

dt 


u - u x p; 


dq> 

dt 


Al 


3) Ihe equation of motion: 




where 



dl ., 

w - — - I', 

dt t' \ 


v being related 

to z by the followit 

:g equation: 


bI K \ 

» ‘ js (z * V’ 

i 
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whence 


All the variables entering into the fundamental system of equations 
can be expressed in terms of the path / and of its derivatives, since 
v - i;, 7 - 


jfm 

si 


expressed as a function of 1^, the pressure is proportional to l^, 

and the derivative — is proportional to Consequently, if the 

dt t 

time t of travel of the projectile through the bore is taken as the 
independent variable, and the path of the projectile / is taken as 
the function to be expanded in a series, it becomes possible to 
employ Taylor’s series for finding the value of the path and 

of its derivatives for the neighboring segment corresponding to the 
time t n+1 - t R + At « t n + h, provided that the values of the path 
l D and of its derivatives for the preceding instant t R are known. 

It is thus possible to find all the elements of burning of the powder 
and of the motion of the projectile during a shot, i.e., z, v, p, 
l, and t. 

Let it be assumed that for a certain instant of time t n the path 
l n and its derivatives with respect to time l^, /jJ, are knownj 

if a sufficiently small Increment of time At - h is assumed and 
consideration is limited to derivatives up to and including the third 
order, then, in accordance with Taylor's formula, we shall have for 
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(108) 


i . * + h !: + “/:♦ — /-. 

n-fl n n 2 n 2»3 n 


Differentiating with respect to t and rejecting the terms 

containing i.e., considering that f"' is constant over the 

n n 

given interval At and equals its mean value, we obtain: 


+ + L /- , 


r + hr' 

n n av. 


j*n + j'*t 

where . ■ - B ills the mean value of the third derivative in the Interval 
2 

under consideration (Fig. 154). From the last equation we obtain the 


1" 2 /’ 

following value for l n+l • - l r 


n+1 - K 'n ' l n • 



Fig. 154 - Diagram for and . 

An has Men shown by P.V. Melentyev, it is sore convenient to 
compute not the derivatives themselves, but, rather, the quantities 

proportional to them, namely: hi’, h 2 l" and ~( ,M . Therefore, by 

h 3 

multiplying by h and by ~, we will obtain the following 
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hi’ . - h|*.+ h 2 l " * I ••• 
ml n n 2 n 


(109) 


D_ / v« - h |" 

2 ( h+l n ‘n+1 


/ 

- **i r - - £ i- 

R ^ n 


( 110 ) 


Comparing these* two equations with the initial equation (108), 
it will be seen that f enters everywhere without a coefficient, /’ 
enters with the coefficient h, /" enters with the coefficient h 2 , 
icient - 

the subsequent computations. 


r; - 


or, multiplying by h z : 


h 2 l" - h 2 S-p. 

t <fm 


( 111 ) 


L : combining the resulting values for the path l and its 
derivatives with the equations of the fundamental system, we obtain 
the totality of formulas necessary for the solution in the following 
form and sequence, which corresponds to the order of their application, 
the constants encountered being designated below as follows: 


Si K 


fw 


to-JL-k s ~ v 

2 ’ 2i«-» v 2 3 ' <fa 4' s 

n P 


F-TS-7327-RE 


CHARLFS A. VEYr.R K C.n i,. i; 



f hi' - hi' + h 2 i" + 1L I'" ; 
n+1 n n 2 n 


(1) 


hi' 


System I 


n+1 h 


z n+l * 7 o 


V n+l * *0 + Vn+l* 


+»+i ‘"Vn + * X *L 


'on' ' tt + h 'i + T h 2 | i' + i^'; ,! 


v 2 

U, _ Jill 

T n+1 2 . . 2 

_ m _ Yr P . ^n+1 “ k 3 V n+l 

p n+l - * *•> 


**n+l + ^n+1 a +n+l + *.n+l 


- 2| n " + i - » 2 ~ - n+ i * A 4 *„i» 


^ I f,f B 1.21 «t k 2 j «i h |««i 

T *n+l h n+1 * h ‘n • 2 ! n 


(II) 


(III) 


(IV) 


(V) 


(VI) 


(VII) 


(VIII) 


The subscript (n+1) designates those values of the derivatives 
at the end of the given interval of tine which, in the process of 
conputation, are transferred fron the column being computed into the 
corresponding rows of the right-hand neighboring column; however, 
the transferred values now bear the index n because they characterize 
the initial value of the given quantity in the next column. 
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In orc!i*r to perform the computation by means of this totality 
of forrulns, the values t>f / and of its derivatives at the start 
of the projectile's notion, i.e., at the instant t - 0, must be 
known. Since at the start of motion the path I and the speed v 
are equal to 7ero, we obtain: 

<f) 0 - 0, ( D 0 - (v) Q - 0, , ( (” ) 0 - k 4 p 0 , h 2 (|") 0 - h 2 k 4 p o , 


where Pq is the forcing pressure usually specified beforehand. As 
regards the third derivative we shall first find an expression 

for it at the present instant in the form of . 

To determine I"’, we differentiate the equation - k 4 p with 
respect to t: 


I, + I 


^ 4 n-i- V<14 •> 


for the start of motion when I - 0, v - 0, p - p n , arid - M^^we 
will therefore obtain: (p’) ft - ° U 


\ 

Po Vo ' 


t'0 1. 


p 0 \ u 

j ~o 


> 7 T V/ 

r”MhV’ 


The quantity |l + j-£- ] “ and this expression is therefore 

+0 
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sometimes given in the following form: 


Consequently, the value of the third derivative for the initial 
instant is likewise known: 


T'o - WK'o - V V* -f? 


and it is possible to begin the successive solution of System (I) 
first for the first column corresponding to the interval of time 
At « h, then for the second column, etc., thus obtaining a successive 
series of values for /, v, z, *V and P aB functions of t. 

The quantity h - At must be so chosen as to obtain 10-15 columns 
for the period of burning of the powder, which will give a 
corresponding number of points for each of the quantities p, v, /, 


Since the time of burning is fundamentally determined by the 
thickness of the powder, the interval of time h - At may be taken 
approximately according to the formula: 


where e^ - the thickness of the powder in millimeters, h being 
rounded off to one or two significant figures (to 5 in the second 
significant figure). For example, if 2e^ - 1.28 mm: 


t:. A. U‘ 'MfMA 
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h - O.Oul (U.jtU) - O.Utwfi-l - O.OUOti or t».Ut'i‘o5 

Since v A and / A ,irc known iit least approxima U-I v in .iJv.iiire, 

it is | ossl ble, alter computinn the 'average line 1 o! motion of the 
2 / 

r,ro iec tile t * -, to take lor the; value of the time step (increment) 

' ‘ t av ‘ v A 

At - h ~ rounded off to two significant fibres. 


Sequenc e of Computation. All the constants are computed first: 



k . ?w . k _ iiL. k 
1 s i„' 2 s ’ 3 


2fw 2 


tf n 


*<», 


’-v ■4 k *T;r« -h ^T;i i( 

M> 


it 


1/ 
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The sequence ol computation is not ntfectcd regard loss ol whether 
the computation is performed for a duress 1 ve or a progressive' 
powder. 

In computing the segment after theytiecomj os i t ion of the 
progressive powder,it is necessary to substitute for the usual 
formula the following previously derived formula: 

'V- 'Vg + H (z - l)/“l + + ^2^2 (z “ 

/ 

P 

where z varies from 1 to 1 + ——, and and X_ are characteristics 

& & 

of the powder form after decomposition. 

A form for conducting such computations is presented on pages 
691-692. 
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Table 7-c - Computation Form for U.e Solution of Problems ir 
ballistics by Fxjansion in 'laylor’s Series. 


Subscri | t n : 


put at ion Formulas 


a • "'n • h "';; 4 V» 


v - h -in 

N n + 1 h 


0 1.1322 0.3K56 


0.0958 0.1928 0.3546 


0.0364 0.0606 0.1012 


| 0.1322 0.3856 0.8414 


! 0.0308 0.0899 1 

0.0297 0.0297 io.0297 


r n+l * * Xz Ll 


■ 2,850,000 
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12 *3 V n-t 1 


0.0605 0.1196 


0.0641 0.1268 


0.0002 0.0009 


10.0639 0.1259 


0.0002 0.0016 


13 *n+l “ k 3 v Ll 0.0637 0.1243 

14 k 2 ( +nfl “ k 3 V n+l ) " 1S1500 354300 

A n+1 

691 4 V «. 









-C - Computation Form for the Solution of Problems in Internal 
Ballistics by Expansion in 'laylor’s Series. 

Subscri pt n‘ _ 01 23 




Column No. 

1 

2 

3 

Renarks 



lime t R+1 - (n + i )h 

0.0008 

0.0016 

0.0024 



, 1 

h, i 

0 

1.1322 

0.3856 

From line 4 of precod ; 
column 

i 

i 2 

h 2 l» 

hf,.., 

0.0958 

0. 1928 

0.3546 

From line 25 of preceding 
column. Into lii.e 26 of 
the given column. 


3 

0.0364 

0.0606 

0.1012 

From line 28 ol preceding 



2 n 




column. 







Into line 27 of the given 
column. 


4 


0.1322 

0.3856 

0.8414 

Into lii.e 1 of next column 







Into line 16 of next 
column. 


5 

v . h, ^i 

y n + l - 5 - 

165.3 

482 





k l v n+l 

0.0308 

0.0899 





+ 






7 

z o 

0.0297 

0.0297 

0.0297 

In all columns. 


8 

z n+l 

0.0605 

0.1196 




9 

* 7 n + l 

0.0641 

0.1268 





+ 






, 10 

* Xz n+i\ 

0.0002 

0,0009 




11 

1 

+ n+l 1 

O'. 0639 

O. 1259 




12 

k 3 v Ll 1 

0.0002 

0.0016 




13 

*n+l " k 3 v n+i 

0.0637 

0. 1243 




14 

V^n+l " k 3 V h+l ) " 

181500 

354300 





“ A n+1 \ 






691 




l)h O.OOOK 0.0016 0.0024 

0 0.0600 0.3088 

0 0.1322 0.3856 

0.0479 0.0964 0.1773 

0.0121 0.0202 0.0337 


1 

^n+1 

-1 

0.0600 

0.3088 

0.9054 

+ 

3.016 

3.016 

3.016 

+ 'a 

3.076 

3.325 


1 

v J68 

0.134 


- 'n+l ' 'a - 

3.008 

3. 191 


Vl 




- —rc*-*kg, cm 5 

B n + 1 

604 

1100 


h2 'n.l 

0.1928 

0.3546 

- 

h 2 l" 

n 

-0.0958 

-0.1928 

-0.3546 

b? r 

2 n 

-0.0364 

-0.0606 

-0.1012 

b 3 .... 

2 ‘n-H 

0.0606 

0.1012 












n 

0 

1 

2 a 



| Column No. 

1 

2 

3 

Remarks 


Time t , 
n+1 

- (n + l)h 

0.0008 

0.0016 

0.0024 

• 

15 


In 

0 

0.0600 

0.3088 

From lino 19 of 
preceding column 

16 

- 

h 'i 

0 

0.1322 

0.3856 

From line 4 of preceding 
column. 

17 

+ 


0.0479 

0.0964 

0.1773 

(of line 2 of the 

2 

given column) 

18 


"T *2*0 

0.0121 

0.0202 

0.0337 

~ (of line 3 of the given 
3 

column) 



0.3088 0.9054 

3.016 3.016 


Into line 15 of next 
column 


In all columns 




604 

1100 

0.1928 

0.3546 

-0.0958 

-0.1928 

-0.0364 

-0.0606 

0.0606 

0.1012 



_ 

Into lines 2 and 26 of 


next column 

-0.3546 

From line 25 of preceding 


column 

-0.1012 

From line 28 of preceding 


column 


Into lines 3 and 27 of 


next column 















0.0121 


0.0202 


0.0227 


J - ±\ 1.065 I 2 1 j 


0.0600 0.2088 0.9054 


! 3.016 2.016 3.016 


22 


n + 1 n+1 ‘a 


24 P n+1 - ~ ~kg/cro“ 


k 4 h p n+l 


3.076 3.325 


3.008 3.191 


604 1100 


0.1928 0.3546 


-0.0958 -0.1928 -0.3546 


-0.0364 -0.0606 -0.1012 


0.0606 0.1012 


Formulas for the second period: 


v ' w 1 ■ \7~rr 


(1 - k,v ); p 




V - V A P “ P. 
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0.0121 0.0202 



Into lines 2 and 26 of 
next column 


-0.0958 -0.1928 -0.3546 From line 25 of preceding 

column 


M1.0364 -0.0606 -0.1012 From line 28 of precedi:.g 

column 


0.0606 0.1012 


Into lines 3 and 27 of 
next column 


\ Formulas for the second period: 

when l-L v « v. p - p . 


are determined from firsr period when ^ - 1. 
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The extreme left column contains the "computation lormulas" 
and constants of System (I); in the next column to the right these 
lormulas are broken down into individual operations, which are 
followed in the computations. 

To start with, the first column (No. 1) corresponding to the 

time interval 0 to h is filled in first. In this column thc^ubscript 

n relates to the start of the interval, and the subscript n +1 relates 

to its end; for this column n - 0 and n +1 - 1 . 

For the next (second) column, n • 1, n+1 - 2, etc. For the 

first column, computation of the constants gives us at n - 0 I - 0 

n 

(the path at the start of the motion), which we write on line 15; 

h l' - 0 (the velocity at the start of the motion) is written on 
n 2 

lines 1 and 16; h I " - h 2 1" - 0.0958 is written on lines 2 and 26; 

3 n 0 

is written on lines 3 and 27. Line 17 is filled with —(h 2 ("), 
20 i h 3 2 ® 

and line 18 with )• Thus, all the quantities with the subscript 

n • 0 are inserted in the first column. We now subject them,to the 

necessary operations. The sum of the first th'ree rows gives the 

fourth hi' . - hi' which is immediately transferred to lines 1 
n+i 0+1 

and 16 of the neighboring column wherein, provided with the subscript 
n, it characterizes the value of this quantity at the start of the 
next interval; we then determine 7 n+ j» 4 'n+l* and A n+ 1 * 

adding the four rows from line 15 through line 18, we obtain in line 


of time - and this quantity, provided with the subscript n, is 

transferred to line 15 of the neighboring column. After determining 

p +1 - and multiplying it by k.h 2 , we obtain h 2 /" , which we 

"n +1 n+1 

write on line 25 of the first column and on lines 26 and 2 of the 
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next column, where this quantity acquires the subscript n, as does 

1 2 s 
also “(h i" ), which is written.in line 17 of the next column. 

2 n+1 

After performing the operations indicated in the form with lines 

3 

25, 26, and 27, we obtain in line 28 of the first column 

1 Ut 1 

which we immediately transfer to lines 3 and 27 of the next column, 

3 

l v while ) is written in line 18 of the same column. 

• 3 2 n+1 

Thus, all operations with the quantities- bearing the subscript 


t 


V- 


t 


, (■ 


n in the second column are already prepared, and the second column 
is then treated In the sane manner as was the first. 

Constants such as z . in line 7 and in line 20 are inserted 
0 A 

in the series of columns in advance. 

By applying the same rules, to the neighboring second column, 
we shall gradually, step by step, obtain values for v, z, l, 
and p as functions of t - (n+l)h, and this is continued to the end 
of burning or to the instant of emergence of the projectile from 
the bore, it being necessary to use ^- 1 after the end of burning. 

In performing the computations it is necessary to exercise 
extreme care not to commit any errors, because an error in one of 
the preceding columns will distort the results obtained in the 
succeeding columns. 

It is best to follow up the computation of the data in each 
column by plotting them on graph paper as a function of t.. In so 
doing an error in the given column will cause a deviation from the 
regular disposition of the points derived from the preceding columns 
and such an error can be detected and corrected. 

As a criterion of accuracy, it is also useful to plot on the 
h 3 

diagram the third derivative (or 0^ in the last row), which 
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should lirst increase, then pass through a maximum, then become 
zero (p’ - 0) at the instant P m Is attained,and thereupon acquire- 
a negative value, fJuctuating s1ight1y in either direction. 

h 3 

The instant of tiro*' cut oi f on the diagram at p* - 0 or —l"’ - 0 

/ ‘ 2 
corresponds to the iifstant of maximum pressure, and all elements 

for it are best taken from the diagram 

The time t^ corresponding to the end of burning of the powder 
is determined from the diagram on the basis of the 4», t curve at 
4*- 1; thereupon, the elements corresponding to the end of burning 
of the powder for this time are found by interpolation. If, without 
changing the segments At - h, the second period is computed as a 
direct continuation of the first, assuming 4* - 1 and t^*» 1^ throughout, 
the third derivative P" usually begins to fluctuate, sometimes 
entering the region of positive values, which contradicts the 
physical nature of the process of pressure change. 

For this reason, once the elements of the end of burning 
t K , v R , /£, and have been obtained from the computation of the 
first period, the procedure is continued by adopting the same step 
At - h with tjj as the starting point by first computing the values 
of the path and of its derivatives for the start of the second period 
in accordance with the following formulas: 

'(0) - V hl ’(0) - h V h2, '(0, ■ h2 <f» p K - h Vic ! 



t? S. (1 4 
2 <lj + 1 ) 
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whereupon they are written in the initial column for compute ling 
the data of the second period and subjected to the sane operations 
as in the first period, with the sole exceptions that 'l'" 1 is 
assumed in line 11 and that lines 6-10 are omitted. 

The computation is continued in this manner until ^ /^. 

If the remaining elements v A and p A are obtained 

automatically in the same column for the subscript n+1; if ^n+1 >/ A ' 
the computations in this column must be carried as far as line 24, 
with lines 6-10 omitted, whereupon the value of is used to obtain 
t A by interpolation in the last column for the purpose of subsequently 
obtaining the elements p A and v A . 

Instead of expanding in a series after obtaining the elements 
corresponding to the end of burning of the powder, it is possible to 
compute v A and p A by means of the usual second-period formulas, but 
without determining the time t. 

The solution by expansion in a series is applicable to both 
the geometric and the physical law of burning of powder. In the 
latter case: 


* 



from which we have the following expression for I: 



and the correlation between vy and z is replaced by the graphical 
dependence of I upon which is found from the bomb test. 
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The equation for is replaced by the; following expression: 


<P I , 0 * k 2 r o j f- p O 
+0 


and: 


.3 

b v , 
2 n 




where is the value of the experimental function T corresponding 
to the quantity ^ 

In the case of ballooning powders, is greater than in the 
case of the geometric law of burning, and therefore the values of 
both the first derivative and of p itself will increase more 
rapidly, and the maximum pressure will occur earlier. If the 
propellant force of the powder is the same, the maximum pressure will 
be greater in the case of the physical law of burning with ballooning 
than in the case of the geometric law of burning. 
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